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Phase transitions in 1/4-filled quasi-one-dimensional molecular conductors are studied theo- 
retically on the basis of extended Hubbard chains including electron-lattice interactions coupled 
by interchain Coulomb repulsion. We apply the numerical quantum transfer-matrix method to 
an effective one- dimensional model, treating the interchain term within mean-field approxima- 
tion. Finite-temperature properties are investigated for the charge ordering, the "dimer Mott" 
transition (bond dimerization) , and the spin-Peierls transition (bond tetramerization) . A co- 
existent state of charge order and bond dimerization exhibiting dielectricity is predicted in a 
certain parameter range, even when intrinsic dimerization is absent. 
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insulator transition, charge ordering, Mott transition, spin-Peierls transition, bond dimeriza- 
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Quasi-one-dimensional (QID) molecular conductors 
have been intensively studied and new phenomena are 
continuously found, even in well-studied families such as 
DCNQI2X and TM2X with monovalent closed-sheU X 
units. ^ Here, DCNQI is an abbreviation of i?ii?2-DCNQI 
(Ri, R2: substituents) and TM stands for either TMTSF 
or TMTTF. Cationic X+ in DCNQI2X and anionic X' 
in TM2X result in 1/4-filled 7r-bands of DCNQI and TM 
molecules, respectively, in terms of electrons and holes. ^ 
Despite being isostructural within each family, both with 
columns of stacked DCNQI/TM molecules, and having 
similar noninteracting band structures with the same fill- 
ing, they exhibit a rich variety of electronic and electron- 
lattice coupled phases. These facts clearly point to the 
interplay of different interactions, which is of keen inter- 
est in the research of strongly correlated systems. 

The DCNQI salts, where the one-dimensional (ID) 
stacking of DCNQI is uniform at high temperatures (T) , 
undergo various phase transitions to symmetry-broken 
states by cooling. For example, in DI-DCNQl2Ag, charge 
ordering (CO) phase transition at T = 220 K is followed 
by antiferromagnetic phase transition at 5 K.'^"^ On the 
other hand, in DMe-DCNQl2Ag, bond dimerization at 
100 K and bond tetramerization at 80 K, the latter inter- 
preted as the spin-Peierls (SP) transition, are suggested 
by X-ray diffraction measurements.^'^ Varying i?2, 
and X+ , and applying pressure to each compound sensi- 
tively change such properties. Most systems are consid- 
ered to fall roughly into either of the two cases above, 
but still many works need to be carried out to obtain a 
general view. 

Analogous symmetry-broken states are found in 
TMTTF2Ar as well, which are situated in a more cor- 
related regime than their selenide analogs TMTSF2Ar 
because of their smaller bandwidths.* In this family, an 
intrinsic dimerization exists even from high T due to 
the crystal structure, in contrast with the DCNQI com- 
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pounds. CO phase transition is found in many members 
of TMTTF2X,9'i° while at low T either SP or anti- 
ferromagnetic transition takes place. By applying pres- 
sure, the CO phase is suppressed while the interplay be- 
tween CO and magnetism shows a complicated behav- 
ior, ^-'^ which requires revisions of the generic phase dia- 
gram of TM2X.^ 

Numerous theoretical works have been aimed at these 
QID systems.i2 xhe 1/4-filled ID extended Hubbard 
model, considering not only the on-site U but also the 
nearest-neighbor Coulomb interaction V, and its dimer- 
ized version have been studied, as well as electron-lattice 
coupled models. However, still important factors are 
lacking to understand the experimental systems espe- 
cially at finite T, since most of the works are intended to 
clarify the ground-state properties of such ID models. 

Several have actually studied finite-T properties 
of these models. Mean-field studies have been per- 
formed, ^^"^^ but such an approach fails, in general, to 
reproduce paramagnetic insulating phases which are ob- 
served in the above compounds at intermediate T. More- 
over, when fluctuation effects are properly taken into ac- 
count, purely electronic ID models do not show phase 
transitions unless T = 0.^^"^^ On the other hand, finite-T 
phase transitions to electron-lattice coupled phases have 
been investigated by Sugiura et al.,^^^'^^ while recently 
the relevance of interchain Coulomb repulsion has been 
pointed out by Yoshioka et al^^ to reproduce the CO 
transition at finite T. Both of these works were performed 
by means of bosonization, and hence, properties below 
transition temperatures are difficult to describe within 
this method. 

In this letter, toward providing a unified view of the 
various phase transitions at finite T observed in QID 
molecular conductors, a comprehensive study of a model 
including both electron-lattice couplings as well as in- 
terchain Coulomb interactions is presented. We use a 
numerical method which takes into account ID fiuctu- 
ations, playing crucial roles in the phase transitions, and 
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which enables us to study finite-T properties both above 
and below transition temperatures. 

Our Hamiltonian is given hy H ^ X^j I^ehm ~^ + 
+ 'H±, where the first three terms represent each 
chain with index j (we omit it in r.h.s. below). The first 
term reads 



'••EHM 



+ h.C.) 



i.s 



which is the ID extended Hubbard model where standard 
notations are used.^^ We allow intrinsic dimerization 
in the transfer integrals to account for TM2X. The next 
two terms, 



2 



KP2 

4 



(2) 



(3) 



are the Peierls-type and Holstein-type electron-lattice in- 
teractions with coupling constants {gp, -ftTp, Ky'2\ and 
{gHj ^h}j respectively. The lattice distortions Ui and Vi 
are treated as classical variables. Ui represents a change 
in the bond length between sites i and z + 1, measured 
from its equilibrium value in the high-T limit. Vi is the 
lattice degree of freedom affecting the on-site potential 
energy that the electrons experience, which originates 
from, e.g., the e-mv coupling^'^ and/or the coupling with 
the anions.^** Note that we have included the fourth-order 
term with respect to Ui in eq. (2) to avoid physically un- 
reasonable values of gpu^ < —1. We set t to unity and 
also gp = 5h = 1 so that they are incorporated in the 
definitions of Ui and Vi. The chemical potential /i is ad- 
justed to fix the electron density at 1/4. 

The term expresses the interchain Coulomb repul- 
sion between nearest- neighbor chains, Vj_, and we do not 
consider interchain hoppings here. We treat in the 
mean- field approximation following Yoshioka et al.^^ By 
considering possible antiphase charge disproportionation 
between all nearest-neighbor pairs of chains, naturally 
favored by V±, it can be expressed as 



= -zVi_ E 



(4) 



where z is the number of nearest-neighbor chains. Then 
the problem is now reduced to an effective ID problem 
with (rii) determined self-consistently, considering uni- 
form Ui along the interchain direction. 

Thermodynamic properties of the effective ID model 
are calculated by means of the numerical quantum 
transfer-matrix method. ^^"^^ Using the largest eigen- 
value of the transfer-matrix, we have performed calcula- 
tions in the thermodynamic limit at each Suzuki- Trotter 
number m, where more precise results can be obtained 
as m increases. The results, presented below, show that 
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Fig. 1. (Color online) (a) Charge ordering phase transition for 
U = 10, y = 4, and zV± = 1, without the Peierls coupling term, 
and (b) dimer Mott phase transition for U = 10, V = zV±^ = 0, 
and Kp = Kp2 = 0.5. From top to bottom, T-dependences 
of order parameters nco and u^, charge susceptibility Xc and 
magnetic susceptibility Xs are shown. Cases without the order- 
ings arc also shown by open symbols for comparison. The shaded 
areas (in all figures) represent the temperature region where our 
numerical method is less reliable. 

systematic m dependences and physically reasonable re- 
sults are obtained except for low T ( ^ 0.2 t, as indicated 
by shaded areas in the figures shown below). The lattice 
distortions Vi and Ui are self-consistently determined to 
minimize the free energy. As to Vi, the self-consistent 
relation is given by Vi = (gn/Kfi) {ui). Eliminating Vi 
in eq. (3), we find that this term just renormalizes the 
value of zV± in eq. (4) as zV± zV± + g^/KB_; the Hol- 
stein coupling enhances charge disproportionation. In the 
following, we write the renormalized value as zV±_ and 
omit the Holstein coupling term. On the other hand, the 
self-consistent relation for Ui is a cubic equation of Ui 
involving the expectation values, {c] ^Cj^_^_i ^ + h.c). In 
the present calculations, we assume two-fold or four-fold 
periodicity along the chain to seek possible symmetry- 
broken states, and obtain self-consistently converged so- 
lutions for {{ni),Ui}. All the results shown below for 
symmetry breakings with a period of two remain sta- 
ble even when the calculations are extended to four-fold 
periodicity, at least for 0.2 t where our numerical 
scheme is reliable. 

First, let us discuss the results for Sd = 0. Figure 1 
shows typical behaviors in (a) CO and (b) "dimer Mott" 
(DM) phase transitions, both for U = 10. The CO 
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Fig. 2. (Color online) (a) Phase diagram in the region of com- 
petition between the charge order (CO) and dimer Mott (DM) 
states, for U = 10, V = i, and Kp = Kp2 = 0.5. DM+CO de- 
notes the coexisting state of the two orders; schematic drawings 
of the three symmetry-broken states are shown in the figure. The 
thick solid lines are phase boundaries at m — > oo expected from 
finite-m data, (b) T-dependences of CO and DM order parame- 
ters for zV± = 0.7 and 1.1. 

and DM states are stabilized when {V,Vj^} and the 
Peierls couphng dominate, respectively. In Fig. 1(a), be- 
low the CO transition temperature, Tco(— 1-04), the 
charge density at each site becomes disproportionated 
as (ui) — 1/2 + (— l)*nco, where nco is the order pa- 
rameter. At TcOj the charge susceptibility Xc {—dn/d^, 
n: average electron density) shows a kink and decreases 
rapidly below Tco- This shows an opening of a charge 
gap, indicating a metal- insulator transition at Tco- In 
contrast, in the magnetic susceptibility Xs {=dm/dh, m: 
average magnetic moment, h: magnetic field coupled to 
m), we find no noticeable sign at Tco - On the other hand, 
as shown in Fig. 1 (b) , the DM transition is characterized 
by bond dimerization, i.e., Ui — (— l)*Ud, where Ud is the 
order parameter. Finite Ud makes the system effectively 
1/2-filled; therefore, together with U, it spontaneously 
produces a Mott insulator out from the original 1 /4-filled 
band.^* Again, Xc shows a kink at the transition temper- 
ature Tdm(— 0.91) and a decrease below it, suggesting 
a metal-insulator transition. The gradual decrease com- 
pared with the case of the CO transition in Fig. 1(a) 
reflects the slow development of the charge gap in the 
DM transition. 

By tuning the parameters, we can investigate the com- 
petition between the two ordered phases above, i.e., CO 
and DM. This is summarized in Fig. 2(a) as a phase di- 
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Fig. 3. (Color online) Successive dimer Mott and spin-Peierls 
transitions for U = 5, V = V± = 0, and Kp = Kp2 = 0.5. 
(a) T-dependence of lattice dimerization u^, tetramerization ut, 
charge disproportionation nt with a period of four and (b) mag- 
netic susceptibility Xs- The spin-Peierls state is schematically 
drawn in the inset. 

agram varying zV±. In the small- and large-zVx regions, 
DM and CO phases are stabilized, respectively, while in 
between, for zV± ~ 0.4 ^ 1.2, we have found a coexis- 
tence of the two orders (DM-I-CO). As a result, the phase 
diagram shows a tetracritical behavior. Developments of 
the two order parameters, nco and u^, for two values 
of zV± are shown in Fig. 2(b). For zV± = 0.7, the sys- 
tem undergoes successive phase transitions by lowering 
T, as uniform state DM DM-hCC^*^ On the other 
hand, in the region near the tetracritical point such as 
at zV± = 1.1, a further phase transition at low T is ob- 
served, as DM+CO — » CO. We note that the coexisting 
DM-I-CO phase is not found in previous studies on mod- 
els with Sd = 0, and that, moreover, this state exhibits 
dielectricity, analogous to the CO state in electronic ID 
models under 5d ^ 0.^'''^^'^^ 

Next, we demonstrate that our scheme can repro- 
duce transitions to symmetry-broken states with a pe- 
riod of four, which were found in previous studies at 
T — 0.^'^"'^^ Figure 3 shows an example of bond tetramer- 
ization, i.e., SP transition under the DM state ("dimer- 
ization of dimers"),^'^ for U = 5, a, smaller value than 
in the cases above. In Fig. 3(a), one can see that be- 
low Tdm(— 1-02), another phase transition takes place 
at Tsp(~ 0.3) where the bond tetramerization ut ap- 
pears on top of Ud] Ui shows a four-fold periodicity, as 
Ud, ~Ud + Ut, Ud, —Ud — Ut, and so on [sec the inset of 
Fig. 3(b)]. At the same time, a small amount of charge 
disproportionation with a period of four and amplitude 
rit shows up.'^^ Sudden changes of the order parame- 
ters suggest a first-order nature of this SP transition. As 
shown in Fig. 3(b), Xs exhibits a clear drop below Tsp, 
indicating the existence of a spin gap in the SP state. 

Finally, we discuss the results for i5d 7^ 0. In this case, 
bond dimerization exists in the entire range of T, and it 
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Fig. 4. (Color online) (a) Charge ordering and (b) spin-Peierls 
transition for 5d ^ 0. In (a), U = 10, V = 4, zV± = 1, Kp = 1 
(no ii'p2-term), and 5^ = 0.05, and in (b), J/ = 5, V = VI = 0, 
Kp = Kp2 = 0.5, and 5d = 0.01. 

is given by the sum of intrinsic dimerization Sd and self- 
consistently determined Ud when the Peierls couphng is 
taken into account. Nevertheless, CO transition and SP 
transition can occur as in the 6d — case. In Fig. 4(a), 
we observe the CO transition at rco(— 1-03); Ud is al- 
ways finite and shows a cusp at Tco- On the other hand. 
Fig. 4(b) shows an example of the SP transition, which 
is again first-order-like, analogous to the = case in 
Fig. 3. 

Let us discuss our results in relation with experimen- 
tal works. We should point out first that the tempera- 
ture range we have studied is higher than that in the ac- 
tual molecular compounds; experimentally relevant tem- 
peratures are one order of magnitude less than our en- 
ergy unit t.^^ Therefore, the following comparison will 
be qualitative; nevertheless, we believe that our results 
contain the essential physics. The distinct behavior be- 
tween Xc and Xs at the transition temperatures seen in 
Fig. 1 is characteristic of transitions to insulators due 
to strong electronic correlation.^* These are consistent 
with experiments, e.g., in DI-DCNQl2Ag (refs. 3-5) and 
in TMTTF2X compounds exhibiting CO,2'''36,37 ^j^gj.g 
signs of the transition are observed in the charge sec- 
tor such as transport properties while no clear signs are 
found in bulk magnetic measurements. The DM-I-CO co- 
existence obtained in our numerical calculation has not 
yet been found in experiments, and future experimen- 
tal works on the precise nature of each DCNQI com- 
pound are awaited to explore it. The SP transition un- 
der the DM state in Figs. 3 and 4(b) accounts for DMe- 
DCNQIzAg,^'^ and for TMTTFzPFg under pressure," 
respectively. The first-order-like behavior suggested in 
our results is not reported in the above compounds, but 
in a related material EDO-TTF2PF6, a clear first-order 
phase transition is observed. 

In summary, we have theoretically investigated finite- 
temperature phase transitions in quasi-one-dimensional 
molecular conductors by means of the quantum transfer- 
matrix method accompanied by self-consistent determi- 
nation of lattice distortions and interchain mean fields. 
We have clarified behaviors of various phase transitions 
relevant to DCNQI and TMTTF salts, such as the charge 
ordering, the bond dimerization interpreted as a dimer 



Mott transition, and the spin-Peierls bond tetrameriza- 
tion. A coexistence of charge order and bond dimeriza- 
tion exhibiting dielectricity, even when intrinsic dimer- 
ization is absent, is found for the first time. Although 
the temperature range we have studied here is high com- 
pared with that in experiments, the present results tak- 
ing ID fluctuations into account appropriately capture 
the physics of phase transition phenomena in molecular 
conductors. 
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